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The analysis of the recent experimental data on charged-current neutrino-nucleus scattering cross 
sections measured at MiniBooNE requires relativistic theoretical descriptions also accounting for the 
role of final state interactions. In this work we evaluate inclusive quasielastic differential neutrino 
cross sections within the framework of the relativistic impulse approximation. Results based on the 
relativistic mean field potential are compared with the ones corresponding to the relativistic Green 
function approach. An analysis of scaling and superscaling properties provided by both models is 
also presented. 

PACS numbers: 25.30.Pt; 25.30.Fj; 13.15.-|-g; 24.10.Jv 



O 



> 
m 

p 

rn 
O 



I. INTRODUCTION 

From the analysis of the world (e, e') data at quasielas- 
tic (QE) kinematics, i.e., with non nucleonic degrees of 
freedom playing a minor role, it has clearly emerged the 
validity of scaling arguments for values of the transferred 
momentum q high enough [1-4] . Furthermore, the study 
of the separated longitudinal/transverse (L/T) data has 
shown that whereas the L response superscales, scaling 
violations are visible in the T channel due to effects be- 
yond the impulse approximation (lA): correlations, me- 
son exchange currents (MEC), etc. These results consti- 
tute strong constraints for any theoretical model applied 
to electron scattering reactions; not only scaling and su- 
perscaling behavior should be fulfilled, i.e., independence 
on q and on the nuclear target, but also the model should 
be capable to reproduce the specific shape of the exper- 
imental superscaling longitudinal function that presents 
a significant tail extended to large values of the trans- 
ferred energy, that is, positive values of the superscaling 
variable ip [4]. 

A large variety of theoretical models used in the lit- 
erature satisfy scaling arguments. This is the case of 
the relativistic Fermi Gas (RFG) model that superscales 
by construction [2-4]. Other approaches based on non- 
relativistic reductions also fulfill scaling behavior even 
in presence of final state interactions (FSI) and nucleon 
correlations (see Refs. [5-9]). However, most of them 
lead to scaling functions that do not get asymmetry, be- 
ing in clear disagreement with data. Our past studies 
have shown that there are two basic aspects in the de- 
scription of electron scattering which are necessary to 
compare superscaling properties of different models with 
experimental data: relativity and FSL On the one hand, 
scaling comes out only when the momentum transfer is 
high enough. This means that, not only relativistic kine- 
matics should be considered, but also the nuclear dy- 



namics and the description of the operators should in- 
volve relativity. On the other hand, the importance of 
FSI has also been clearly stated in the case of exclu- 
sive (e, e'N) processes where the use of complex optical 
potentials is required [10-18]. However, even being re- 
stricted to the QE kinematical domain, the analysis of 
inclusive reactions, contrary to exclusive ones, needs all 
inelastic channels in the final state to be retained, i.e., 
the flux should be conserved. Therefore, the distorted 
wave impulse approximation based on the use of a com- 
plex potential should be dismissed for (e, e') processes. 

Previous studies [19-22] have clearly illustrated the key 
role played by FSI to reproduce the specific asymmet- 
ric shape shown by the experimental scaling function. 
Specifically, description of the final nucleon state using 
the same relativistic mean field potential considered in 
describing the initial state, denoted simply as RMF, leads 
to a superscaling function in accordance with data. Re- 
cently [22], we have extended this analysis by introduc- 
ing the relativistic Green function (RGF) technique de- 
veloped by the Pavia group [23-26]. This is an alter- 
native, relativistic, approach to FSI in (e,e') reactions. 
Both descriptions, RMF and RGF, contrary to previous 
models based on non-relativistic or semirelativistic re- 
ductions [5] , provide a significant asymmetry in the scal- 
ing function following the general behavior of data. It is 
also important to point out that asymmetry in the scal- 
ing function can be also obtained within the framework 
of a semirelativistic (SR) model, provided that FSI are 
described by the Dirac equation-based (DEB) potential 
derived from the RMF [9]. 

High-quality predictions for neutrino-nucleus cross sec- 
tions are needed for use in on-going experimental stud- 
ies of neutrino oscillations at GeV energies. Recently, 
muon charged-current neutrino-nucleus double differen- 
tial cross sections have been measured for the first time 
by the MiniBooNE collaboration [27]. The experimen- 
tal conditions have made it possible to disentangle data 
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for which no pions in the final state are detected. This 
corresponds to the QE regime where only nucleons are 
the relevant degrees of freedom, and is simply denoted 
as charged-current quasielastic (CCQE) neutrino-nucleus 
scattering. The analysis of CCQE data shows that the 
REG model understimates the total cross section, unless 
the axial mass Ma is significantly enlarged with respect 
to the accepted world average value. The kinematics in- 
volved in the MiniBooNe experiment: neutrino energy 
flux ranging from 375 MeV up to ~ 3 GeV, requires the 
use of a relativistic description of the process. Moreover, 
a proper analysis of data requires also to have good con- 
trol on nuclear effects [5, 19, 20, 24, 28-43]. Notice that 
any reliable calculation for neutrino scattering should 
first be tested against electron scattering in the same 
kinematical conditions. In this sense, while the REG is 
capable of getting the basic size and shape of the in- 
clusive cross section, it can hardly account for important 
details of the response. Thus, the need for a modification 
of the axial mass within the REG could be a way to ef- 
fectively incorporate nuclear effects into the calculation. 
As mentioned, relativity is expected to play an impor- 
tant role at MiniBoone conditions. Not only relativistic 
kinematics, already included in the REG model, but also 
dynamical relativistic effects, not present in REG, can be 
important. In this work, the two QE-based approaches 
selected, RME and RGE, incorporate relativistic dynam- 
ics by means of strong scalar and vector potentials in the 
Dirac equation employed to describe nucleons in nuclei. 

As shown in Refs. [19-22], the RME and RGE ap- 
proaches describe successfully the behavior of electron 
scattering data and their related scaling and superscaling 
functions. Therefore, in this work we extend our previ- 
ous analysis for (e, e') reactions to charged-current (CC) 
neutrino-nucleus scattering. The differences observed be- 
tween the predictions of the two approaches can be help- 
ful for understanding nuclear effects, particularly ESI, 
which may play a crucial role in the analysis of neutrino- 
nucleus scattering data and its influence in studies of neu- 
trino oscillations at intermediate to high energies. 

The paper is organized as follows. In Section II we 
summarize the basic formalism involved in the descrip- 
tion of CC neutrino-nucleus scattering. We introduce 
the main expressions and present a brief discussion on 
the two approaches considered in the description of fi- 
nal state interactions, namely, the relativistic mean field 
and the relativistic Green function technique. Scaling 
and superscaling arguments for neutrino-nucleus scatter- 
ing processes are also briefly mentioned. In Section III 
we show and discuss the results obtained with the two 
models considered. After checking the reliability of both 
calculations in the relativistic plane wave approach and 
also making use of the real part of the relativistic opti- 
cal potential, we describe ESI within the RME and RGE 
models, and present results for the differential cross sec- 
tions and scaling functions. Einally, in Section IV we 
summarize our conclusions. 



II. NEUTRINO-NUCLEUS SCATTERING 

The single-nucleon cross section for an inclusive re- 
action where an incident neutrino (antineutrino), with 
four- momentum fc'^ = (£'y,fc), is absorbed by a nu- 
cleus and only the outgoing muon, with four-momentum 
k'^ = {e\k') and mass m ~ 105.7 MeV is detected, is 
given by [5, 6, 44, 45] 
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G ~ 1.166 X 10^5 GeV"2 jg ^he Eermi constant, dc is 
the Cabibbo angle with cos^ (i^c) — 0.975, and i? is the 
generalized scattering angle, that is defined by 



tan^ (i9/2) = 
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The four- momentum transfer is = /c^ — k'^ = (w, q), 
with ~ Lo^ — . The factor comes from the 
contraction between the Icpton and the hadron tensor 
and contains all the relevant nuclear structure informa- 
tion, can be decomposed into a charge-charge (CC), 
charge- longitudinal {CL), longitudinal- longitudinal {LL) 
and two transverse {T,T') responses, i.e.. 



■^+{-} — VccRcc + '^VclRcl + VllRll + VtRt 



2xVt'Rt' , 
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where x ~ +1(^1) ii^ the case of neutrino (antineutrino) 
scattering. The coefficients V are given by 

Vcc = , Vcl = - + ^ ts^n" {d/2) , 

= (^)'+(l + |^ + ^^p)^^tanM^/2) , 

V ( 7 + ^ ) tan' ('^/2) , 
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The weak response functions are given by 
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i{W ~W^')/2 , (7) 



in terms of the components of the hadron tensor, that is 
given by suitable bilinear products of the transition ma- 
trix elements of the nuclear many-body charged-current 
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operator between the initial state | ^fo) of the nucleus, 
of energy Eq, and the final states | ^'/), of energy Ef, 
both eigenstates of the {A + l)-body Hamiltonian i7, as 

X {^„\rHq)\^s) 6{E„+uj-Ef),{%) 

involving an average over the initial states and a sum over 
the undetected final states. The sum runs over the scat- 
tering states corresponding to all of the allowed asymp- 
totic configurations and includes possible discrete states. 

In the QE region the scattering is described in the 
one-boson exchange approximation and in the relativistic 
impulse approximation (RIA) by the sum of incoherent 
processes involving only one nucleon scattering, i.e., the 
scattering occurs with only one nucleon, which is sub- 
sequently emitted, while the remaining nucleons of the 
target behave as spectators. Within the RIA the nu- 
clear current operator is assumed to be the sum of single- 
nucleon currents j^, corresponding to the weak CC op- 
erator 

- G^(g2)7V + ^^p(Q')QV]^*, (9) 



where are isospin operators, k is the anomalous part 
of the magnetic moment, and a^^^ = {i/2) [7^,7"]. 
and are the isovector Dirac and Pauli nucleon form 
factors, which are related to the corresponding electro- 
magnetic form factor by the conservation of the vector 
current. Ga and Fp are the axial and induced pseu- 
doscalar form factors which are usually parametrized as 



Ga 
Fp 



9A 



(1 + QVMl) 

2mNGA 
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where gA — 1.267, is the pion mass, ttin the nucleon 
mass, and Ma — 1.03 GeV is the axial mass [46-48]. 

Within the RIA framework and under the assumption 
of a shell-model description for nuclear structure, the 
components of the hadron tensor are obtained from the 
sum, over all the single-particle (s.p.) shell-model states, 
of the squared absolute value of the transition matrix 
elements of the single-nucleon current 



{x^E\E)\riq)\^n) 



(12) 



where Xe ' (F) is the scattering state of the emitted nu- 
cleon and the overlap ipn between the ground state of the 
target | ^0) and the final state | n) of the residual nucleus 
is a s.p. shell-model state. In this work the bound nu- 
cleon states ipn are taken as self-consistent Dirac-Hartree 
solutions derived within a RMF approach using a La- 
grangian containing cr, w and p mesons [49-53]. 



Different prescriptions are used to calculate the rela- 
tivistic s.p. scattering wave functions. In the relativistic 
plane wave impulse approximation (RPWIA), FSI be- 
tween the outgoing nucleon and the residual nucleus are 
neglected. In the relativistic distorted wave impulse ap- 
proximation (RDWIA) FSI effects arc accounted for by 
solving the Dirac equation with strong relativistic scalar 
and vector optical potentials. This approach has been 
very successfuU in describing exclusive {e,e'p) data [11- 
15, 54, 55], but it would be inconsistent for the inclusive 
scattering, where all the inelastic channels should be re- 
tained and the total flux, although redistributed among 
all possible channels due to FSI, must be conserved. In 
the RDWIA (with complex potentials) the flux is not 
conserved and the inclusive (e, e') cross section is under- 
estimated [23, 56-58]. A simple way of estimating the 
inclusive response within the RIA is to use purely real 
potentials. In a first approach, the imaginary part of 
the phenomenological relativistic energy-dependent op- 
tical potentials [59] is set to zero, thus retaining in the 
calculations only the real part. In a second approach, 
the scattering states are described by using the same real 
scalar and vector RMF potentials considered in the de- 
scription of the initial bound state (/?„. We refer to these 
two different FSI descriptions as real relativistic optical 
potential (rROP) and RMF calculations, respectively. 

We note that the rROP conserves the flux and thus 
it is inconsistent with the exclusive process, where the 
final state includes only the elastic rcscattcring of the 
knocked-out nucleon and thus a complex optical poten- 
tial adjusted to elastic proton scattering data must be 
used. Moreover, the use of a real optical potential is 
unsatisfactory from a theoretical point of view, since it 
is an energy- dependent potential, reflecting the differ- 
ent contribution of open inelastic channels for each en- 
ergy. Dispersion relation then dictates that the poten- 
tial should have a nonzero imaginary term [60]. On the 
other hand, the RMF model is based on the use of the 
same strong energy-independent real potential for both 
bound and scattering states, so that it fulfills the disper- 
sion relation [60] and, further, it maintains the continuity 
equation. 

In a different approach. Green's function techniques 
may be exploited to derive the inclusive response. This 
formalism allows to reconstruct the flux into non-elastic 
channels in the case of inclusive scattering, starting from 
the complex (scalar and vector) optical potential which 
describes elastic scattering data. The components of the 
hadron tensor are written in terms of the s.p. optical 
model Green's function. This is the result of suitable 
approximations, such as the assumption of a one-body 
current and subtler approximations related to the im- 
pulse approximation. The explicit calculation of the s.p. 
Green's function is avoided by using its spectral repre- 
sentation, which is based on a biorthogonal expansion in 
terms of a non Hermitian optical potential H and of its 
Hermitian conjugate "H^. Calculations require matrix el- 
ements of the same type as the RDWIA ones in Eq. (12), 
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but involve eigenfunctions of both Ti and Ti^ where the 
different sign of the imaginary part gives absorption in 
one case and gain of flux in the other case. Thus, in 
the sum over n the total flux is redistributed and con- 
served. This relativistic Green's function model allows 
for a consistent treatment of FSI in the exclusive and 
in the inclusive scattering and gives a good description 
of (e,e') data [22, 23]. Detailed discussions of the RGF 
model can be found in Refs. [22-26]. 



A. Scaling at the quasielastic peak 

Scaling arguments applied to neutrino reactions were 
presented for the first time in Ref. [5], where a phe- 
nomcnological "SupcrScaling" approach, denoted as 
SuSA, was proposed based on the assumed universality of 
the scaling function for electromagnetic and weak interac- 
tions. Hence, the phenomenological scaling function ex- 
tracted from the analysis of electron scattering data was 
used in order to make "model-independent" predictions 
for neutrino-nucleus cross sections. The universal scaling 
behaviour is highly related to the single-boson exchange 
approximation and this justifies that both neutrino and 
electron scattering in QE conditions exhibit to a large 
extent similar scaling. A different way to analyze the 
problem consists in making use of a specific model that 
works properly in the description of electron scattering 
data, and apply it to neutrino reactions. In this case, 
although cross sections are "model-dependent" , theoret- 
ical scaling functions corresponding to neutrinos can be 
derived and compared directly with the theoretical func- 
tions for electrons (obtained in the same model) and also 
with the experimental scaling function. Thus, the uni- 
versality property can be tested within the framework of 
a specific model. 

In this work we apply scaling arguments to the RMF 
and RGF models. Results are shown in the next sec- 
tion. The usual procedure to get the scaling function has 
been considered, i.e., dividing the inclusive cross sections 
(I) by the appropriate single- nucleon charged-current i/N 
elastic cross section weighted by the corresponding pro- 
ton/neutron number involved in the process. The scaling 
function /{"tp') obtained with this procedure depends on 
the dimensionless scaling variable ip'{q, w) extracted from 
the RFG analysis that incorporates the typical momen- 
tum scale for the selected nucleus [4, 20]. The explicit 
expressions for the CC singlc-nuclcon responses have al- 
ready been presented in Refs. [5, 6]. In these single- 
nucleon responses, effects coming from the motion of the 
nucleons in the nucleus have been considered only up to 
first order in the Fermi momentum rjp = kp/m^. The 
use of the full results for the non-Pauli blocked regime, as 
shown in [6] , leads to very minor differences without mod- 
ifying the discussion and general conclusions presented in 
the work. 




400 600 800 400 600 800 




800 1000 1200 1400 800 1000 1200 1400 
[MeV] [MeV] 

FIG. 1: Differential cross section of the ^^C{v^, fi~) reaction 
as a function of the kinetic energy of the outgoing muon r,i 
for two values of the incident neutrino energy, i.e., = 1000 
and 2000 MeV, and fixed muon scattering angle = 45 deg, 
calculated by the Pavia (solid lines) and the Madrid-Sevilla 
(dashed lines) groups with RPWIA [(a) and (c)] and rROP 
[(b) and (d)]. 



III. RESULTS AND DISCUSSION 

In this section the numerical results of the different rel- 
ativistic models developed by the Pavia and the Madrid- 
Sevilla groups to describe FSI in the inclusive quasielas- 
tic CC neutrino-nucleus scattering are considered. As a 
first step results obtained by the two groups in the RP- 
WIA and rROP approaches arc compared to check the 
consistency of the numerical programs when calculations 
are carried out under the same conditions. Then the 
results corresponding to the RMF model developed by 
the Madrid-Sevilla group and the RGF model developed 
by the Pavia group are compared. The comparison is 
performed for the ^^6(1^^, /i~) cross section and scaling 
function calculated in different kinematics which are rep- 
resentative of the kincmatical range where QE scattering 
is expected to give the main contribution to the inclu- 
sive neutrino-nucleus scattering. Pion-production events 
may play a significant role and increase their contribu- 
tion with larger values of the energy transfer. However, 
its contribution for the CCQE MiniBone data is expected 
to be very minor: no pions in the final state are detected. 
Antineutrino-induced scattering can also be calculated, 
but it is not considered in this paper. All the calcu- 
lations are performed using the same relativistic initial 
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FIG. 2: Diflterential cross section of the ^^C{u^,^ ) reaction 
for three values of the incident neutrino energy, i.e., = 1000 
(a), 1500 (b), and 2000 MeV (c), and = 45 deg. The sohd 
and long dot-dashed lines are the RGF results calculated with 
the two different optical potentials EDADl and EDAI-12C. 
The dashed lines are the results of the RMF model. 



FIG. 3: Differential cross section of the ^^C{v^,fi ) reaction 
for El, = 1000 (a) and 2000 MeV (b) and scattering angle i?,^ 
= 45 deg. The solid, long dot-dashed, and dashed lines are 
the same as in Fig. 2. The dot-dashed and dotted lines are 
the rROP and the RPWIA results, respectively. 



state wave functions, that are taken as Dirac-Hartree so- 
lutions of a relativistic Lagrangian written in the context 
of a relativistic mean field theory with the NLSH param- 
eterization [51-53]. 

A. Differential cross section 

The ^^C {vfj,, /i^) cross section calculated by the Pavia 
and the Madrid-Sevilla groups in RPWIA and in rROP 
for two values of the incident neutrino energy, i.e., E^, ~ 
1000 and 2000 MeV, and a fixed muon scattering angle 
1?^ = 45 degrees are compared in Fig. 1. Almost identical 
results are obtained in RPWIA. In rROP the two results 
are very similar, up to a few percent. A similar agreement 
between the RPWIA and rROP results of the Pavia and 
Madrid-Sevilla groups was already found in Rcf. [22] for 
the inclusive QE electron scattering. The comparison 
in Fig. 1 is an important and necessary benchmark of 
our independent computer programs, which allows us to 
estimate the numerical uncertainties and gives enough 
confidence on the reliability of both calculations. 

Cross sections evaluated with the RMF and RGF mod- 
els for the kinematics with E„ = 1000, 1500, and 2000 
MeV and = 45 degrees arc presented in Fig. 2. In 
the case of the RGF approach two different parameter- 
izations for the relativistic optical potential have been 



used: the energy-dependent and A-dependent EDADl 
and the energy-dependent but A-independcnt EDAI-12C 
complex phenomenological potentials of [59], which are 
fitted to proton elastic scattering data on several nuclei 
in an energy range up to 1040 McV. In the figures the re- 
sults obtained with EDADl and EDAI-12C are denoted 
as GFl and GF-edai, respectively. We note that all the 
RGF results presented here contain the contribution of 
both terms of the hadron tensor in Eq. (25) of Ref. [24], 
while only the contribution of the first term was included 
in the numerical results presented in Ref. [24]. The sec- 
ond term is entirely due to the imaginary part of the op- 
tical potential and vanishes if only the real part is con- 
sidered. The calculation of this term is a complicated 
and time consuming task that requires integration over 
all the eigenfunctions of the continuum spectrum of the 
optical potential. Indeed the first term gives in general 
the main contribution. The contribution of the second 
term depends on kinematics and becomes more impor- 
tant for higher values of the energy transfer, i.e., higher 
energies for the optical potential. Although small and 
even negligible in many situations, the second term can 
be very important in particular kinematic conditions and 
in general cannot be neglected. 

To make the comparison of the different models clearer, 
we compare in Fig. 3 the RMF and RGF results directly 
with the RPWIA and rROP ones. Wc note that in these 
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FIG. 4: Differential cross section of the ^^C{u^,^ ) reaction 
for E,. = 1000 MeV and q = 500 MeV/c (a) and 1000 MeV/c 
(b). Line convention as in Fig. 3. 



kinematics, that are similar to those actually explored 
at present neutrino experiment facilities, the momentum 
transfer q is not fixed and its value around the peak of 
the cross section is usually large, i.e., q ~ 700 MeV/c for 
= 1000 MeV, g « 1100 McV/c for = 1500 MeV, 
and q « 1400 MeV/c for E^ = 2000 MeV. The shape of 
the RMF cross section shows an asymmetry, with a long 
tail extending towards lower values of the energy of the 
final muon, i.e., higher values of w, which is due to the 
strong energy-independent scalar and vector potentials 
present in the RMF approach. The asymmetry increases 
with larger incident neutrino energies. In the case of the 
RGF cross sections the asymmetry towards higher w and 
E^ is less significant but still visible. Almost no asym- 
metry is found for the RPWIA and rROP cross sections. 
We note that in the present calculations, that are per- 
formed in kinematics different from those considered in 
Ref. [24] , significant differences are obtained between the 
rROP and RGF results. These differences stress the im- 
portant role played by the imaginary part of the optical 
potential, which in the RGF model affects both terms of 
the hadron tensor. The GFl and GF-cdai cross sections 
have somewhat different shapes, that are particularly vis- 
ible for low u! at E^ = 1000 MeV and for higher ut at 
E^ = 2000 MeV. These differences are essentially due to 
the imaginary part of the ROP, that is sensitive to the 
particular parametrization considered. The real terms of 
the phcnomcnological ROP are very similar for different 
parametrizations and give very similar results. 



FIG. 5: L, T, and T' components of the cross section of the 
^^C{iy^,fi') reaction for = 1000 MeV and q = 500 MeV/c 
(panels (a), (c), and (d) in the left column) and 1000 MeV/c 
(panels (b), (d), and (e) in the right column). Line convention 
as in Fig. 2. 



The differences between the results of the GFl and 
GF-edai cross sections depend on the energy and mo- 
mentum transfer, and are directly linked to the specific 
structure of the energy-dependent relativistic optical po- 
tentials adopted in the RGF model. The differences be- 
tween the RMF and RGF results increase with the energy 
and momentum transfer. The larger differences seen for 
the largest value of E^, not only between the RMF and 
RGF models, but also between the two RGF results, is 
simply an indication of the difference in the ingredients 
of these calculations. 

To perform a more direct comparison between our 
neutrino- and electron-scattering calculations of Ref. [22] , 
we present in Fig. 4 our cross sections calculated for a 
fixed value of the incident neutrino energy, E^, = 1000 
MeV, and two values of the momentum transfer, i.e., q = 
500 and 1000 MeV/c. Also in these kinematics the shape 
of the RMF cross section shows an asymmetry with a 
tail extending towards higher values of w. An asymmet- 
ric shape towards higher uj is shown also by the RGF 
cross sections, while no visible asymmetry is given by 
the RPWIA and rROP results. Also in these kinematics 
significant differences are obtained between the RGF and 
rROP cross sections. 

The behavior of the RMF and RGF results as a func- 
tion of q and to is connected to the different relativistic 
potentials used in the models. The RMF model uses the 
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FIG. 6: Scaling function of tlie ^^0(1^;^,/^") reaction in the 
same kinematics as in Fig. 2 compared with the averaged 
experimental scaling function. Line convention as in Fig. 2. 



FIG. 7: Scaling function of the ^^C{vi_i,fi ) reaction for inci- 
dent neutrino energy = 1000 MeV and q = 500 MeV/c (a) 
and 1000 MeV/c (b). Line convention as in Fig. 4. 



effective mean field, that reproduces the saturation be- 
havior of nuclear matter and the properties of the ground 
state of nuclei. On the other hand, the RGF model uses 
the phenomenological optical potentials fitted to elastic 
proton- nucleus scattering. The lost of elastic flux into 
inelastic channels caused by the imaginary term of these 
potentials is recovered for inclusive scattering in the RGF 
making use of dispersion relations. As already shown for 
(e, e') reactions [22], the RGF yields a larger cross section 
than the RMF. The latter can be considered as an esti- 
mation of the the pure nucleonic contribution. On the 
contrary, the RGF may account for in a phenomenolog- 
ical way additional effects due to non nucleonic degrees 
of freedom, more important with increasing transferred 
momentum. Notice that both the RGF and the RMF 
yield predictions within few percent for low-g (panels (a) 
of Figs. 2 and 3), while as q goes up the RGF yields 
increasingly larger cross sections than RMF. This may 
reflect the influence of the pionic degrees of freedom. 

The results obtained with GFl and GF-cdai are con- 
sistent with the general behavior of the phenomenolog- 
ical energy-dependent relativistic optical potentials and 
the differences are basically due to their different imagi- 
nary part, that includes the overall effect of the inelastic 
channels and is not univocally determined by the elastic 
phenomenology. 

The results displayed in Fig. 4 present some differences 
with respect to the corresponding cross sections of the in- 
clusive electron scattering shown in Ref. [22], which are 



calculated with the same models for FSI and at the same 
values of the momentum transfer. In both cases the dif- 
ferences between the results of the different models are 
generally larger for increasing value of the momentum 
transfer. For neutrino scattering, however, such a behav- 
ior is less evident and clear. In particular, the GFl cross 
section at q = 1000 Mev/c docs not show the strong en- 
hancement in the region of the maximum that was found 
for the (e, e') calculations of Ref. [22], where the GFl re- 
sult was even larger than the RPWIA one. In the case of 
neutrino scattering the RGF results in the region of the 
maximum are generally larger than the RMF ones, but 
smaller than the RPWIA cross sections. 

In spite of many similarities, inclusive electron scatter- 
ing and CC neutrino-nucleus scattering are two different 
processes and caution should be drawn on their compar- 
ison. We note that the cross sections shown in Fig. 4 
are calculated with same incident lepton energy and the 
same momentum transfer as in the (e, e') calculations of 
Ref. [22]. Nevertheless the muon mass gives in the case of 
CC neutrino scattering a different kinematics, with dif- 
ferent values of the energy transfer and, as a consequence, 
of the energies of the outoing nucleon. This means that 
in the RGF model the optical potential is calculated for 
electron and neutrino scattering at different energies. To 
clarify this point and reproduce the kinematics of elec- 
tron scattering, we have performed some calculations for 
the {v, reaction with vanishing muon mass. The main 
difference with respect to the calculations shown in Fig. 4 
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FIG. 8: Analysis of first kind scaling, for E„ = 1000 

MeV and q = 500 MeV/c (solid lines) and 1000 MeV/c (dotted 
lines) with the GFl (a), GF-edai (b), and RMF (c) models 
using the same results already displayed in Fig. 7. 



is a shift of the cross section, by about 100 MeV, towards 
higher values of , without any significant change in the 
shape or in the strength. 

The most important difference between electron- and 
neutrino-scattering at the Feynman amplitude level is 
the exchanged vector boson: a virtual photon probing 
the electromagnetic current in electron scattering, and 
a probing the weak current in CC neutrino-nucleus 
scattering. Thus, whereas for the (e, e') process all nu- 
cleons (protons and neutrons) in the nucleus contribute, 
in the case of CC neutrino (antineutrino) scattering only 
the neutrons (protons) in the target contribute to the in- 
clusive cross section. Moreover, whereas the Rl and Rt 
responses in CC neutrino induced scattering are given by 
the sum of pure vector- vector (VV) and axial-axial (AA) 
terms (see, e.g. Refs. [5, 6]), the i?T' response comes 
from the interference between the vector (V) and axial 
(A) terms in the neutrino CC operator. In the electron 
scattering case, where both isoscalar and isovector form 
factors enter the responses, the cross section is the sum 
of the longitudinal and transverse responses. 

The different currents and their possible interplay with 
the other ingredients of the models do not allow an easy 
comparison between the results of electron and neutrino 
scattering. The numerical differences between the RGF 
results for electron and neutrino scattering can mainly be 
ascribed to the combined effects of the weak current, in 
particular, its axial term, and the imaginary part of the 



relativistic optical potential. Wc have checked that these 
effects can explain the fact that in neutrino scattering 
the GFl result does not give the strong enhancement in 
the region of the maximum of the cross section that was 
found for the (e, e') calculations of Ref. [22]. 

As a further example, we compare in Fig. 5 the RMF 
and RGF results for the longitudinal, transverse, and in- 
terference axial- vector contributions to the cross section. 
The sum of these three terms gives the corresponding 
cross sections in Fig. 4. Here, by longitudinal contribu- 
tion to the cross section we mean the sum of the Rll-, 
RcL, and Rcc responses, defined in Eq. 7, each one mul- 
tiplied by the appropriate kinematical factors of Eq. 5. 
In general, the RMF and RGF results as a function of 
q and lo show the same behavior as the cross sections in 
Fig. 4 and, again, are related to the different relativistic 
potentials used in the models. The T and T' components 
are similar for all cases: at g = 500 MeV/c both RGF 
results are larger than the RMF ones in the peak region, 
while at g = 1000 MeV/c the RMF and GF-edai resuhs 
are very similar and lower than GFl. The longitudinal 
cross section has a different behavior: at g = 500 MeV/c 
the GF-edai result is close to RMF and lower than GFl, 
while aX q = 1000 MeV/c it overlaps GFl and they both 
are larger than RMF. At g = 500 MeV/c all the GF 
and RMF calculations give the T cross section approx- 
imately two times larger than the T' one. In addition, 
the longitudinal contribution, although smaller, cannot 
be neglected; in this case the relative difference between 
the GFl and GF-edai longitudinal components is large 
and has visible effects in the cross sections of Fig. 4. At 
q = 1000 MeV/c the T and T' components have similar 
strengths, whereas the longitudinal response gives only a 
very small contribution to the cross section. In the elec- 
tron scattering case, &t q = 500 MeV/c the longitudinal 
and transverse components have similar strength whereas 
at g = 1000 MeV/c the transverse response is much more 
important (see, e.g., Refs. [22, 23]). 



B. Scaling functions 

The effects already discussed for the differential cross 
sections are also present in the scaling functions. Here we 
compare results for the scaling function extracted from 
the differential cross section /(■)/'') using the same descrip- 
tions for the final state interactions already considered for 
the cross sections. The transverse responses Rt and Rt' 
arc the leading terms for the total cross section whereas 
the longitudinal contribution R^ is usually small; the 
corresponding scaling function friip') and fT'i'4'') are 
very similar, up to a few percent, and, moreover, they 
are similar to the scaling function fiip') from the cross 
section, i.e., scaling of zero kind is verified. Thus, we do 
not show /t(V'') and fr'iip') and present our results only 
for f{ip'). The longitudinal contribution leads to a scal- 
ing function /l ('(/'') which may depart significantly from 
f{ifi')- However, one should be cautious because the Ion- 
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gitudinal contribution to inclusive neutrino-nucleus cross 
sections can be almost negligible in some kinematics com- 
pared with the transverse T and T' ones. 

As a first step, we have compared the scaling function 
fiip') obtained by the Pavia and Madrid-Sevilla groups 
in RPWIA and rROP for three values of the neutrino 
energy, = 1000, 1500, and 2000 McV, and z?^ = 45 
degrees. Our results arc almost coincident and, in ad- 
dition to the cross section results in Fig. 1, confirm the 
consistency of the numerical codes when calculations are 
performed under the same conditions. We do not show 
these results for simplicity. 

In Figs. 6 and 7 we compare the scaling function 
evaluated with different models and for the same kine- 
matical conditions as in Figs. 2 and 4, respectively, with 
the averaged QE phenomenological longitudinal scaling 
function extracted from the analysis of (e, e') data [2-4]. 
As already shown in previous works [19-21], the RMF 
model produces an asymmetric shape, with a long tail 
in the region with ?/)' > 0, that follows closely the be- 
havior of the phenomenological function. It has been 
noticed that the RMF approach is capable of explaining 
the asymmetric behavior of the phenomenolgical scaling 
function within the framework of the relativistic impulse 
approximation taking advantage of its strong relativistic 
scalar and vector potentials. The results with the RGF 
model in Fig. 6 are similar to those obtained with RMF 
at El, = 1000 McV, while visible discrepancies appear 
at E^ = 1500 and 2000 McV, and in Fig. 7. As a gen- 
eral remark, these results for the scaling functions follow 
similar trends to those already applied to the behavior 
of the cross sections in Figs. 2 and 4, i.e., the differences 
between the RMF and RGF, and between the GFl and 
GF-edai calculations, increase with larger energy and mo- 
mentum transfer. 

The asymmetric shape with a tail in the region of posi- 
tive V'' is obtained in both RMF and RGF models, which 
involve descriptions of FSI cither with a strong energy- 
independent real potential or with a complex energy- 
dependent optical potential, respectively. The scaling 
functions corresponding to RPWIA and rROP, which are 
also presented in Fig. 7 for two values of the momentum 
transfer, do not show any significant asymmetric tail for 
ip' > 0. 

It has been already discussed [19-21] that the RMF 
model produces a scaling function f{ip') that is in ac- 
cordance with the longitudinal scaling function extracted 
from the (e, e') reaction and hence with the electron scat- 
tering longitudinal data. This outcome reinforces the 
validity of the general assumption [20] to predict CC 
neutrino-nucleus cross sections from the phenomenologi- 
cal scaling function, and it is a noticeable fact that this 
is contrary to what one would expect, since neutrino re- 
actions are totally dominated by the transverse Rt and 
Rt' responses. 

In the case of the RGF model we apparently do not ob- 
tain a similar good agreement with the longitudinal scal- 
ing function extracted from the (e, e') results of 



Rcf. [22]. This is clearly related to the fact that, at higher 
q values, there is not such a strong enhancement of the 
maximum strength, particularly for GFl as in Ref. [22]. 
However, taking advantage of this, the RGF model pro- 
duces a scaling function in reasonable agreement with the 
averaged experimental scaling function. 

An analysis of the scaling of first- kind, i.e., indepen- 
dence of the momentum transfer, is illustrated in Fig. 8. 
The results are the same already shown in Fig. 7, but are 
presented in a different way. Each panel corresponds to a 
specific description of FSI (GFl, GF-edai, and RMF) and 
includes the results obtained at g = 500 and 1000 MeV/c. 
The experimental (e, e') data are compatible with a mild 
violation of the first-kind scaling, particularly in the posi- 
tive ?/;'-region. In Refs. [19, 20] the scaling functions eval- 
uated with the RPWIA and rROP models were shown to 
depend very mildly on the transferred momentum in the 
whole, positive and negative, ip' region. In the case of 
the RMF approach, there is a slight shift in the region 
-0' < (slighter than for the longitudinal scaling function 
of electron scattering) , whereas the model breaks scaling 
approximately at 30% level when -0' > 0. Similar re- 
sults arc obtained with the RGF models, where a small 
shift in the region of negative ^' also occurs, and scaling 
is broken for ip' > Q. This scaling violation for > 
is slightly larger with GFl than GF-edai. However, the 
quality of the comparison between the experimental QE 
scaling function and the two RGF results is similar. Fi- 
nally, the stronger scaling breakdown shown by RGF, 
compared to RMF, may refiect the effective presence of 
non-nuclcnic contributions in the case of RGF. 



IV. SUMMARY AND CONCLUSIONS 

This work emerges from the collaboration between the 
Pavia and Madrid-Sevilla groups, and it extends the 
study presented in Ref. [22] where we focused on inclusive 
electron scattering reactions. Here we consider the case of 
charged-current quasielastic neutrino-nucleus scattering 
processes. These studies are of great interest in connec- 
tion with the recent data measured at MiniBooNE [27], 
and their "possible" impact on neutrino oscillations. 

The kinematics of the MiniBooNE experiment involves 
neutrino energies up to ~ 3 GeV. This implies that rela- 
tivity is an essential ingredient; not only relativistic kine- 
matics should be considered, but also nuclear dynamics 
and current operators should be described within a rela- 
tivistic framework. The analysis performed by our groups 
involves a relativistic description of the reaction mecha- 
nism with strong scalar and vector potentials, thus pro- 
viding a proper treatment of the neutrino scattering pro- 
cess of interest. However, before entering into a detailed 
analysis of the MiniBooNE experiment and its compari- 
son with our theoretical results, that will be considered in 
a subsequent work, in this paper we restrict ourselves to 
a systematic comparison between the two approaches for 
different kinematics. As already mentioned in Rcf. [22], 
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the description of the final state interactions constitutes 
a basic aspect in the analysis of the theoretical results 
and their comparison with data. 

Following Ref. [22], we check the consistency of the 
numerical calculations by comparing results in the plane 
wave limit (RPWIA) and also making use of the real 
part of the rclativistic optical potentials (rROP). Almost 
identical results are obtained in RPWIA and very close 
in the case of rROP. This outcome is consistent with our 
previous analysis in electron scattering, and it reinforces 
the reliability of the calculations. Next we apply the two 
models considered in Ref. [22] to account for FSI, that 
is, the RMF, based on the use of the same relativistic 
mean field potential for the bound and ejected nucleon 
wave functions, and the RGF approach, which treats FSI 
consistently in the inclusive and exclusive reactions. In 
the latter we have considered two difl'erent optical poten- 
tials, EDADl and EDAI-12C, which have been fitted to 
proton elastic scattering data. 

Contrary to RPWIA and rROP approaches, the differ- 
ential cross sections obtained within RMF and RGF show 
a significant asymmetry with a tail extended to large val- 
ues of the energy transfer (small muon kinetic energies). 
The amount of this asymmetry in the RMF approach in- 
creases with larger neutrino energies. On the contrary, its 
contribution is less significant in the RGF model, show- 
ing also important differences attached to the two optical 
potentials considered. These discrepancies are linked to 
the imaginary energy-dependent term in the ROP that is 
very sensitive to the particular parameterization consid- 
ered. It seems that optical potentials fitted only to elas- 
tic scattering data are not well constrained when applied, 
within the RGF approach, to describe inclusive reactions. 
In addition to these effects introduced by the particular 
ROP considered, also significant differences emerge from 
the comparison between RMF and RGF, being larger for 
increasing transferred momentum. Although this general 
behavior was already present in the case of (e, e') pro- 
cesses [22], the results for neutrino reactions show some 



distinct features to be discussed. In particular, the GFl 
{v,ii~) cross section does not present the strong enhance- 
ment in the region of the maximum that was previously 
found for (e, e') aX q = 1000 MeV/c. This result has im- 
plications for the scaling function. As shown in Ref. [21], 
the function obtained for CCQE neutrino reactions 

within the RMF model is in accordance with from 
(e, e'). This is connected to the isospin content in the nu- 
cleon form factors. On the contrary, such an agreement 
does not occur for RGF. This would indicate that the 
universality property of the superscaling function is not 
entirely fulfilled within the RGF model. 

Summarizing, we present a systematic analysis of 
cross sections and scaling properties for CCQE neutrino- 
nucleus scattering processes. We consider two different 
relativistic descriptions of the final state interactions. Al- 
though both approaches respect scaling and superscaling 
behavior, some differences emerge, particularly concern- 
ing the comparison between the results obtained for elec- 
trons and neutrinos. The present results should be cau- 
tiously viewed due to the very different ingredients con- 
sidered by RMF and RGF models; however, these may 
help in disentangling different physics aspects involved in 
the processes, with a special mention to the analysis of 
the MiniBooNE data and FSI effects. Work along this 
line is presently in progress. 
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